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Color superconducting matter in a magnetic field 
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We investigate the effect of a magnetic field on cold dense three-flavor quark matter using an 
effective model with four-Fermi interactions with electric and color neutrality taken into account. 
The gap parameters Ai, A2, and A3 representing respectively the predominant pairing between 
down and strange (d-s) quarks, strange and up (s-u) quarks, and up and down (u-d) quarks, show 
the de Haas- van Alphen effect, i.e. oscillatory behavior as a function of the modified magnetic 
field B that can penetrate the color superconducting medium. Without applying electric and color 
neutrality we find A2 — A3 S> Ai for 2eB > fi^, where e is the modified electromagnetic coupling 
constant and is one third of the baryon chemical potential. Because the average Fermi surface 
for each pairing is affected by taking into account neutrality, the gap structure changes drastically 
in this case; we find Ai ^ A2 — A3 for 2eB > /i^. We point out that the magnetic fields as strong 
as presumably existing inside magnetars might induce significant deviations from the gap structure 
Ai ~ A2 — A3 at zero magnetic field. 

PACS numbers: 12.38.Aw, 12.38.-t, 24.85.-|-p, 26.60.-l-c 



By analyzing Quantum Chromodynamics (QCD) it 
has been established that at zero temperature and high 
enough baryon densities color superconducting (CSC) 
matter should be formed jlj. Unfortunately no experi- 
mental evidence for color superconductivity is yet avail- 
able. In CSC matter Cooper pairs of quarks are created 
due to an attractive interaction between quarks on oppo- 
site sides of the Fermi surface. The (almost) sole place 
where one might be able to find color superconductivity 
in nature would be the central part of neutron stars. To 
this aim one has to clarify the properties of CSC matter 
under the physical conditions maintained inside neutron 
stars 

The neutron star density is at most p ^ lOpo in 
the core, where po is the normal nuclear density ^ 
0.17nucleon/fm'^. This density roughly corresponds to 
a quark chemical potential (i.e. one third of the baryon 
chemical potential) /Zq ^ 500 MeV as deduced from p ^ 
3/Xq/7r^. At this intermediate density one cannot neglect 
the role of the strange quark mass Ms = 100^200 MeV. 
The strange quark mass induces a "pressure" to tear the 
Cooper pairs apart, i.e., a Fermi surface mismatch of size 
Mg /2^q between u/d quarks and s quarks will be formed. 
The pairing pattern is quite complicated in the density 
region where /2pq is comparable to the gap energy A 
which are both of order tens MeV around the region of 
our interest where pq ^ 500 MeV. 

The Fermi surface mismatch is further caused by the 
requirement of neutrality which is broken by Mg 7^ in 
three-flavor quark matter. The system should be elec- 
tric neutral to avoid divergent field energies faster than 
the volume, otherwise the system is not stable thermo- 
dynamically. Regarding color neutrality the constraint 
is more stringent, that is, the whole system must be a 
color-singlet. To consider the phase structure, however, 



it is adequate to impose global color neutrality as well as 
global electric neutrality. In the effective model we will 
use in this Letter this can be achieved by introducing 
the electric and color chemical potentials, pe, MSi and 
/ig corresponding the the negative electric charge ma- 
trix Qe = —Q = diag(— |, i, i) in flavor space and two 
diagonal color charge matrices T3 = diag(i,— i,0) and 
Tg = diag(i,i,— |) in color space These chemical 
potentials mimic gauge field dynamics. 

In neutron stars another source for the "pressure" on 
Cooper pairs is a magnetic field. In regular neutron 
stars the magnetic field strength on the surface is of or- 
der B = 10^ ^ 10^^ G and it reaches values as large as 
B - 10^5 G (in which case eS = 1 - 10 MeV^) in a spe- 
cial kind of neutron star called magnetar [^). Actually 
the virial theorem Q enables us to deduce B < 10^* G 
in the interior of the neutron star. In compact stars 
that are self-bound rather than gravitationally bound, 
the maximum magnetic field could be even larger such 
that eB « p'^. Clearly, the effect of these magnetic fields 
are not to be neglected at all as compared to A and 

Mf/2pq. 

Much work has been done to investigate the effect of 
the magnetic field on nuclear matter (see Ref. Q and 
references therein) and also some on norm al q uark mat- 
ter 0,11] but only hmited results [1, [l^jElj 121 are avail- 
able on CSC matter in a penetrating magnetic field. 

In particular in Ref. [13] the analytical solution of the 
gap equation is found only in the limit of a strong mag- 
netic field without taking account of the neutrality condi- 
tions. The new material we shall elucidate in this Letter 
is twofold; first, we solve the gap equations numerically in 
order to get the gap parameters for any value of magnetic 
field within the framework of the Nambu-Jona-Lasinio 
(NJL) model. This model correctly describes qualitative 
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features of high-density QCD [l|, Next, we impose 
the electric and color neutrality conditions on the sys- 
tem which changes the qualitative behavior of the gap 
parameters significantly from the non-neutral case. 

We assume only the predominant pairing in the color 
anti-symmetric channel; 

where a and (3 run from 1 to 3 in color space {r,g,b) 
and i and j run from 1 to 3 in flavor space (u^d,s). The 
gap parameters Ai, A2, and A3 represent the ds-pairing, 
su-pairing, and wd-pairing, respectively. 

This gap pattern breaks the electromagnetic U(l) sym- 
metry. As a result the photon becomes massive, so that 
a pure electromagnetic field cannot penetrate CSC mat- 
ter. However, the pairing pattern is still invariant under 
U(1)q transformations, where Q = Icoior'SJQ — QO^lflavor- 
The corresponding rotated electromagnetic field is = 
cos 9 — G^sin^, which is a combination of the elec- 
tromagnetic field Af^ and a component of the gluonic 
field Gj^. The rotated photon stays massless, hence a 
rotated magnetic field B can penetrate CSC matter Q. 
The coupling constant for A is e = e cos 9 where e and 
g are the electromagnetic and QCD coupling constant. 
Here the mixing angle 9 may depend on the gap struc- 
ture 0, [13] . In the convention that the gauge fields are 
defined with generators normalized as tr[tai6] = 25afc, 
one finds cos9 = g/{^e^ + g'^Y^^ for color-flavor locked 
(CFL) matter (Ai w A2 ~ A3). For the phases with 
only A2 or A3 nonzero (i.e. 2SCsu or 2SC phase) one 

has cos^ — 5/(j^e^ -I- g^)^^^. Interestingly enough, we as a function of = 
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FIG. 1: Gap parameters as a function of eB / fj,q for /iq 
500 MeV without neutrality. 
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For the charged sectors a constant magnetic field re- 
sults in the Landau quantization of the energy dispersion 
relation. For instance, the rs-bu sector of the Nambu- 
Gor'kov propagator contains sixteen dispersion relations 
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Once wc turn on 



find that the mixing angle in the phase with only Ai 
nonzero (i.e. 2SCds phase) is the same as that in the 
CFL phase [H. 

By careful consideration of the boundary layer between 
CSC matter and normal quark matter one can derive a 
relation between the magnitude of the applied external 
magnetic field -Bcxt outside and the rotated magnetic field 
B inside CSC matter. For sharp boundaries (boundary 
smaller than screening length) a small part of the flux is 
expelled and one has B « i?cxt cos 6*, while for smooth 
boundaries B k, iJ^xt hence no flux is expelled at all. 9j. 
Since g 3> e in the region we are interested in, cos 9 k, 1 
which implies that the magnitudes of the magnetic flelds 
outside and inside the CSC core are approximately equal, 
that is, we will implicitly assume eB « eBoxt in our dis- 
cussions. 

While the quark Cooper pairs are neutral with respect 
to the Q charge, some of the individual quarks which 
form a pair are, however, charged under the B field. In 
particular, the three flavors and three colors result in 
nine different quarks from which four are Q charged. We 
summarize in the following table the quark species, Q, 
and the gap parameters involved in the pairing; 



a constant B in the z-direction, rotational symmetry is 
broken and the dispersion relations are modifled by the 
following replacement. 
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where n = 0, 2,3. . . with ±|ei?| depending on the spin 
and Q charge. Then, in order to write down the ther- 
modynamic potential $7 of the NJL model, we can utilize 
the conventional expression (see e.g. Ref. [l6| including 
full account of Ms ^ 0) with the prescription for the 
rs-bu and gu-rd sectors. The associated momentum in- 
tegration is modified as follows 
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We will present the prove of this simple replacement in 
a separate paper [3]. 

We first solve three gap equations without neutrality 
taken into account, 
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at He = M3 = Ms = to check consistency with Ref. [lo| . 

We show our numerical results as a function of a dimen- 
sionless parameter eB/ix^ in Fig. [T] for = 500 MeV. 
We chose the cut-off parameter A = 1000 MeV and the 
four-Fermi coupling constant to yield Aq = 50 MeV for 
Ms ^ B = 0. We made sure that A-dependence is 
tiny once the coupling constant runs as a function of 
A to give a fixed value of Aq. To reduce the cut-off ar- 
tifact, we used a smooth Fermi-Dirac-like form factor, 
^{1 — tanh[{pB — A)/a;]} with a choice uj = 0.1 A in the 
momentum integration. If we adopt a smaller value of uj 
which results in a sharper cut-off scheme, the curves in 
Fig. [1] become less smooth. In particular, for u; — > we 
find tiny spikes in the gap parameters when h? / {2eB) is 
an integer, originating from the vacuum energy contribu- 
tion to the thermodynamic potential. We have checked 
the robustness of the smooth oscillatory shapes seen in 
Fig. [1] by varying the value of a; [31 ■ 

We can see that A2 and A3 are close to each other 
apart from a discrepancy by the strange quark mass 
which is Ms = 100 MeV in our calculation; the effect 
of Ms pushes Ai and A2 down to 47.7 MeV and A3 up 
to 50.8 MeV at B = 0. The gap parameters show oscilla- 
tory behavior as long as 2eB < because the density of 
states increases every time 2eBn approaches /Zq for some 
non-zero n. Thus, the oscillation ceases when the first 
Landau level (n — 1) lies above the Fermi surface, i.e. 
2eB > /iq. In fact, it is manifest in Fig.[T] The behavior 
of the gap parameters is similar to the oscillation of the 
magnetization of a material in an external magnetic field 
which is known as the de Haas- van Alphen effect. 

Our numerical results are qualitatively consistent with 
the analytical evaluation in Ref. [l3| for eB ^ /Xq when 
only the lowest Landau level (LLL) contributes to the gap 
equations. In that case the gap parameters A2 and A3 
increase monotonically as a function of B. The reason 
for this is understood in view of the analytical expres- 
sions given by Eqs. (95) and (96) in the second paper of 
Ref. [1^ ; roughly speaking, as a result of the LLL approx- 
imation, the phase space is enlarged as 2^q — > 11"^ + eB. 
At the same time we see that Ai decreases slightly as 
a function of B because it is only indirectly sensitive to 
J5, which is qualitatively consistent with Eq. (101) in 
Ref. [l3l. 

Next we will impose three neutrality conditions, 
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dfie 9^3 9^8 
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(5) 



We have also added the contribution of a free electron and 
muon gas to f2. Electrons and muons feel the magnetic 
field BcosO with the coupling constant e, that amounts 
to eB cos 9 = eB. Once we take account of neutrality the 
situation drastically changes. Figure [2] shows the gap pa- 
rameters with fie, /i3, and /ig determined self-consistently 
for /iq = 500 MeV. The corresponding chemical poten- 
tials are displayed in Fig. [3] It can be seen in Fig. [3] 
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FIG. 2: Gap parameters as a function of eB/fj,^ with neutral- 
ity, for = 500 MeV. 



that at very large eB, fi^ is larger than the muon mass 
~ 106 MeV. This indicates that in that case a significant 
number of muons will be present in the system. 

In sharp contrast to the non-neutral case, we find that 
Ai grows with increasing eB. The gap parameters A2 
and A3 vanish smoothly at eB / fjL^ ~ 1.09 and 1.47, re- 
spectively. This implies a second order transition from 
the CFL phase to the so-called dSC phase, followed by 
a transition to the so-called 2SCds phase. Note that 
the mixing angle 9 is common in the CFL, dSC, and 
2SCds phases. The behavior of Ai at large eB can be 
accounted for by /i3 and /ig; the Fermi surface average 
P'bd-gs becomes larger, for example in our calculation, 
from 500 MeV at i? = to 625 MeV at eB/^l = 4. This 
results in a larger gap parameter. So the system exhibits 
a phase with only Ai nonzero at large eB/fj,^. This im- 
plies two-flavor color superconducting pairing between d 
and s quarks, hence the name 2SCds. The possibility of 
the 2SCds phase as a ground state for eB/fi^ > 1.47 
is quite interesting, since the 2SCds phase has rarely 
been paid attention to in the QCD phase diagram. [See 
Ref. 
gion. 



17[ for detailed analyses including the 2SCds re- 



The 2SCds phase is similar to the more familiar 2SC 
phase where only A3 takes a finite value. The behavior 
of /i3 and /ig in the 2SCds region might look totally dif- 
ferent from the 2SC phase in which fj,^ = and /ig <C /ie- 
This difference, however, turns out superficial once we re- 
arrange the color-flavor bases properly as fJ-'^ = f^s ~ ^ fJ-s 
and /ig — —5/^3 — /ig — |/ie. Then, as shown by dotted 
curves in Fig. [31 /tg in the 2SCds phase is zero just like /i3 
in the 2SC phase, and /ig stays smaller than /ig just like 



4 



100 



> 

4^ 50 



S 

o 

a 

o 

n 

2-50 



-100 



y V /» .•• •• 




— Me 

- /i3 



■•■ - M3/2 

-■ -/i3/2 - ^8 - 3Me/2 

0.1 



FIG. 3: Chemical potentials as a function of eB/ fi^, for fiq 
500 MeV. 



fig in the 2SC phase. We note that /X3 is not oscillatory 
at all even though fi^ and ns are. We also point out that 
CFL matter remains a Q-msulator, i.e. /ig = [l8| until 
ei?//^2 ^ 0.88. 

In fact, there are two windows in which gapless dis- 
persion relations with non-zero gap parameters appear; 
0.88 < eS/^q < 1.09 for rs-bu pairing with A2 and 
1.25 < eB/y?^ < 1.47 for gu-rd pairing with A3. The 

2SCds phase for eB / > 1.47 is fully gapped. Once 
the system enters the gapless state, A2 and A3 rapidly 
decrease to zero. Because the phase space is enlarged 
by large B for Q-charged quarks and thus their density 
increases, the Fermi surface averages firs-bu and ftgu-rd 
should be located lower than the others to keep neutral- 
ity. It can happen with developing /13 ^ and fi^ 
which alter not only the Fermi surface averages but also 
the Fermi surface mismatches 5^rs-bu and Sjigu-rd- Under 
the constraint of neutrality, hence, A2 and A3 become 
smaller because of the decreasing Fermi surface average 
with increasing _B, and at the same time, the Fermi sur- 
face mismatches for A2 and A3 grow up with /i3 and 
/ig induced by B, and eventually the gapless dispersion 
relations emerge when S^rs-bu > A2 or S^igu-rd > A3. 

The question of whether the 2SCds phase is a real pos- 
sibility of the ground state under a sufficiently strong 
magnetic field or not should be answered by energy com- 
parison with normal quark matter. We have to calculate 
the energy density in CSC matter with eB and that in 
normal quark matter with eB^xt to determine which is 
energetically favored. Unfortunately the energy density 
has a huge oscillation as a function of B if eB ~ . This 
oscillation arises from the vacuum energy contributions 



to the thermodynamic potential. Because normal quark 
matter couples in another way to the magnetic field, the 
oscillation in the energy density of normal quark mat- 
ter is different from that of CSC matter, making energy 
comparison at eB ~ ambiguous. 

We have carefully checked that the cut-off dependence 
is tiny in the gap equations and neutrality conditions as 
we mentioned, and it is natural because the momentum 
integration near the Fermi surface should be dominant. 
Therefore we could as well have taken a larger value of 
A such that the gap parameters remain to have the same 
value and eB is much smaller than A^. In this way we 
found that there is always an energy gain between CSC 
and normal quark matter, in particular it is kept to be 
at least 2 x 10® MeY* around eB/^l - 1. 

However the field energy it takes to expel part of the 
applied magnetic field from CSC matter is in the case 
of a sharp boundary equal to ^(-Boxt ~ cos^flSpxt) ^ 
e^i3gxt/6(7^. The coefficient e^/g^ is of order 0.01, and 
therefore, the energy cost is of order 1 x 10® MeV* around 
eB/fi^ ^ 1, which is comparable or less than the energy 
gain. In the case of a smooth boundary no flux is ex- 
pelled at all, so then CSC matter is always favored. We 
might thus expect to see Ai 3> A2 ~ A3 before reaching 
a possible transition to normal quark matter. Hence if 
ei?cxt ~ Aiq, the 2SCds phase could be a likely candidate 
for the ground state of matter inside magnetars. 

In summary, we considered the effect of a strong mag- 
netic field on neutral CSC quark matter. We found that 
the neutrality conditions significantly change the non- 
neutral results. We pointed out the possibility of the 
2SCds phase in the interior of the magnetar. 
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